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Abstract 

Thermal fluctuations and disorder play an essential role in high-Tc supercon- 
ductors. After reviewing the mean-field phase diagram we describe significant 
modifications that result when the effects of finite temperature and disorder 
are incorporated. Thermal fluctuations cause the melting of the Abrikosov flux 
line lattice into a vortex liquid at high temperatures, and disorder produces 
novel glass phases at low temperatures. We analyze phase transitions between 
these new phases and describe transport properties in various regimes. 
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1 INTRODUCTION 



The microscopic mechanism of superconductivity in the new high-Tc materials has 
remained a mystery despite much vigorous recent activity However significant 
progress in our understanding of the behavior of these materials has been achieved 
through phenomenological approaches. The phenomenological theory of supercon- 
ductors was introduced in the 1950's by Landau and has been successful in describing 
the conventional superconductors [Q. Many of the key properties and phenomena 
in high-Tc materials can also be described entirely with the conventional Ginzburg- 
Landau theory, albeit with unconventional values of parameters. In this new range 
of parameters an understanding of the interplay of strong thermal fiuctuations and 
disorder, both in statics and dynamics, is required. 

Within a mean-field treatment, appropriate to conventional (low-Tc) type II su- 
perconductors, the Ginzburg-Landau theory leads to the celebrated Abrikosov vor- 
tex lattice 1^ and the Meissner superconducting phases, with well-defined continuous 
transitions occurring at the upper and lower critical fields Hc2(T) and Hci(T), respec- 
tively (see Fig.||(a)). However, unlike these conventional type II materials, where in 
the absence of disorder the flux lines always arrange themselves into an ordered, static 
triangular lattice, high-T^ superconductors exhibit much richer vortex-line states. Be- 
cause of the combination of high Tc (~ 100 K), short coherence length ^, large pene- 
tration length A, and large anisotropy, the Abrikosov lattice of flux lines is destabilized 
at flnite temperatures by the thermal fluctuations and disorder. 

The resulting magnetic fleld versus temperature (HT) phase diagram for high-Tc 
materials is drastically modifled by the fluctuations and disorder. In Fig.|l|(b) we 
have schematically illustrated a resulting phase diagram for the high-Tc materials, 
with various new novel phases and regimes that have been proposed. The fluctuation 
and disorder-corrected phase diagram is remarkably rich and is to be contrasted with 
mean-field phase diagram characterizing the conventional superconductors, displayed 
in Fig.|I|(a). 

At low temperatures, the translational order of the Abrikosov lattice is lost due 
to disorder that appears in the form of oxygen vacancies and interstitials, and grain 
and twin boundaries 0]. Vortex glass 0, Bose glass and polymer glass [0 are 
some of the new low-temperature phases that have been proposed to replace the 
traditional Abrikosov fiux lattice. Even in pure crystals thermal fiuctuations have 
been proposed to lead to a novel supersolid vortex phase ||, which unlike the 
3d Abrikosov lattice, is a quasi-2d lattice with freely wandering interstitials. At 
even higher temperatures or higher fields the fiux lattice melts and is replaced by a 
vortex liquid of highly mobile and flexible vortex lines |TD|, |TT[ . Motivated by recent 
experiments, in which a resistive "shoulder" was observed, proposals have been made 
that the vortex liquid phase could possibly be further subdivided into distinct liquid 
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phases, and can undergo transitions between hexatic [O, O, H and isotropic, and 



entangled and disentangled vortex liquids |10 



To study the transport properties in these new novel phases the interaction be- 
tween the transport current, the screening supercurrents and the magnetic field must 
be understood. These interactions lead to a Magnus-Lorentz force exerted on the vor- 
tex line perpendicular to the transport current. A motion of vortices in response to 
this force generates finite voltages and resistivity, and therefore leads to a breakdown 
of true superconductivity. 

The Meissner phase, in which the fiux is completely expelled, is fully supercon- 
ducting up to a critical surface current density jc beyond which thermal fiuctuations 
can produce vortex loops, which expand under the infiuence of the transport current. 
In the ideal situation an unpinned vortex lattice can move as a whole in response to 
a uniform current and therefore is not truly superconducting. However, the vortex 
lattice has no resistive linear response for periodic current patterns with a finite wave- 
length because of the finite shear modulus of the lattice. Therefore the vortex lattice 
exhibits a novel fully nonlocal linear resistivity. Of course in practice the vortex lat- 
tice will have a superconducting linear response even to a uniform current, since the 
sample boundaries and even weak disorder will pin the lattice in place. Newly pro- 
posed vortex glass phases are also believed to be fully superconducting. The vanishing 
of linear resistivity in the Meissner and the putative glassy solid phases (vortex and 
Bose glass) is related to the divergence of energy barriers (with vanishing current) 
separating the ground state and the low-lying excited states. 

In the vortex liquid phase the fiux lines can more easily respond to the transport 
current-induced forces, although interactions with other lines are believed to lead to 
a much higher viscosity than that of the conventional isotropic liquids [0, 0. A 
phenomenological hydrodynamic theory can be used to describe the resulting fiux- 
fiow resistivity which is nonlocal with properties intermediate between that of the 



fully superconducting vortex lattice and the normal phase |T5|, ^ . 

The mean-field phase transitions are also modified by strong thermal fiuctuations. 
In contrast to low-T^ materials, in the Cu02 based superconductors there is no well 
defined thermodynamic transition at the upper critical field, Hc2- The mean-field 
transition is replaced by a smooth crossover (indicated by a dashed line in Fig.[I|(b)) 
from the normal state to the non-superconducting, but more conductive vortex liquid 
regime. Therefore the distinction between the vortex liquid and the normal state is 
only quantitative. As H^2^(T) is passed from above the resistivity decreases, the 
vortices begin to form as the fiuctuations in the superconducting order parameter 
grow to a nonlinear regime. However, because of the motion of vortex lines, which 
destroys phase coherence, the average value of the superconducting order parameter 
?/'(r) is still zero in vortex liquid regime as it is in the normal state. 

The real transition occurs only at lower fields and temperatures, where supercon- 
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ducting vortex solid phases form. For clean samples a first order freezing transition 
transforms the liquid into the vortex lattice. Despite the lack of a full theory of this 
3d melting transition the location of the melting curve Bm{T) can be estimated by 
the Lindemann criterion [|l^, |l^, |l^]. However, in the presence of disorder and strong 
vortex line interactions physical arguments have predicted a transition to the highly 
correlated fully superconducting isotropic and anisotropic (Bose glass) vortex glass 
phases. The vortex glass transition is believed to replace the melting transition of 
pure samples. Although no fundamental analytical description is available to date, 
this continuous glass transition can be analyzed using scaling theory that is in good 
agreement with transport experiments. 

Recently, there have been several experimental observations of the transition into 
a glassy vortex phase that agree with predictions of the scaling theory |2y. 



pT| , p2| . The impressive data collapse over several decades, from which the critical 
exponents characterizing the transition can be determined, provides good evidence 
for the existence of the vortex glass phase. This data and analytical analysis is further 



supported by various numerical studies [^, 0, ^ 



In these lectures we will review some of the highlights of the phenomenology of 
high-Tc superconductors. We begin in Sec. 2 by reviewing the picture that emerges 
from the mean-field analysis, describing the equilibrium and non-equilibrium proper- 
ties of the resulting phases and the mean-field transitions. As we described above, 
effects of thermal fluctuations and disorder are very important and in Sections 3 and 
4 we describe how these effects modify the HT phase diagram. After describing equi- 
librium properties of the new vortex phases, we proceed to analyze their transport 
properties. We apply a scaling theory to predict the form of the resistivity in the 
superconducting phases. In Sec. 5 we describe the nature of the fluctuation-corrected 
phase transition and again use scaling theory to summarize its properties in terms of a 
small set of universal critical exponents, which can be determined through numerical 
and experimental means. In each section we compare theoretical results with exper- 
iments and numerical studies, and for the most part find a coherent and consistent 
picture emerging. 



2 MEAN-FIELD THEORY 

Irrespective of the microscopic mechanism responsible for the attractive interaction 
that leads to the binding of the electrons into Cooper pairs, the Ginzburg-Landau 
theory postulates the existence of an order parameter, a complex scalar field ilj{f) 
that is the condensate "wavefunction" of the bosonic Cooper pairs of electric charge 
2e. The long wavelength properties of the superconductors are then encoded in the 
Ginzburg-Landau free energy functional that is assumed to have a local expansion 
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in powers of iplf) and its spatial gradients. Near the superconducting transition, Tc, 
the order parameter ipl-f) is small, and it is sufficient to retain only the lowest order 
terms in the expansion of the free energy [EH 




The first two terms are the supercurrent kinetic energy. The sum of the third and 
fourth terms is the local pairing energy, approximated at quartic order. The last 
term is the magnetic field energy that couples the external magnetic field to the 
magnetic induction i?(r^ = V x A{f). Physically the parameters and m in the 
above expansion can be identified with the effective mass of the Cooper pair along 
the c-axis (chosen in 2;-direction) and in the a6-plane. With the anisotropy parameter 

7 = ^m/niz ^ 1 we can describe the large material anisotropy of the high-Tc layered 
materials, which are much stronger superconductors in the afe-plane than along the 
c-axis 1^^. The a and (3 are the material parameters that depend on temperature 
T. Generally a(T) ^ a{T — T^^^)/T = t while the other parameters have only 
a smooth dependence on T without sign changes. This dependence and the general 
validity of above description has been verified by Gorkov |28[ through his derivation of 
the Ginzburg-Landau theory from the microscopic Bardeen-Cooper-Schreiffer (BCS) 
theory of conventional superconductors 0. 

The mean-field description which neglects the effects of fluctuations is valid only 
outside the critical region around the normal-superconducting (NS) transition. For 
conventional, low-temperature bulk superconductors in zero magnetic field, this crit- 
ical region is unmeasurably small with fluctuations only becoming important for 
{T —Tc) /Tc < 10~^. On the other hand the combination of high-Tc, strong anisotropy, 
short coherence length ^ and large magnetic penetration length A of high-Tc materials 
expands the critical region by several orders of magnitude. As we will see, thermal 
fluctuations will lead to significant modifications of the mean-field description that 
we first present below. 
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Figure 1: (a) Mean- field phase diagram of type II superconductors, (b) Schematic pic- 
ture of a phase diagram of high-T^ superconductors which includes effects of thermal 
fluctuations and disorder. 

2.1 Mean-field statics 

The mean-field theory description is obtained by simply minimizing Fgl['4^,A] with 
respect to ijj{r) and A{r). The resulting mean- field phase diagram is displayed in 
Fig. 0(a), with three distinct phases. For T > i.e. a{T) > the pairing and 
magnetic terms are minimized hj ip{r) = and 5(r) = H characteristic of the normal 
state. In the Meissner phase the U{1) original symmetry of Fql ~^ ip{f)e^'^) is 

spontaneously "broken" by ip{r) = = ^J\a\/ [3 

{ip*{r)tp{r + R)) ps , for ^ oo , (2) 

which leads to "off- diagonal" long-range order that sets in exponentially over the 
coherence length ^ = Ti/ \j2m\a\ in the afe-plane. Along the z-axis the coherence 
length is much shorter, = 7^, scaled down by the anistropy parameter, which can 
be as small as 1/100 in BSCCO. 

Borrowing the notation from superfiuids we associate the square of the condensate 
amplitude with the density of the superconducting pairs. It plays the role of stiffness 
for the phase 0(r) variations, as we will see below. The coherence length over which 
ipif) varies (due to for example imposed boundary conditions) is set by the compe- 
tition between the kinetic energy and the pairing energy and reflects the size of the 
Cooper pair. 

The magnetic field i?(r) = is expelled everywhere in the sample except for the 
boundary where it leaks in a distance A = ^Jmc'^j3/16n\a\e'^ along the afe-plane. The 
superconducting screening currents along the 2;-direction are less effective in keeping 
the applied field from penetrating along the c-axis, with A2 = A/7. The London 
penetration length A is set by the competition between the magnetic field term and 
the kinetic energy due to the screening currents. 

For vanishing applied field the transition is in the same universality class as that of 
paramagnetic-ferromagnetic transition described by the XY model, with disordered 
paramagnetic state and spin aligned correlated state corresponding to the normal and 
Meissner states, respectively. 

As was predicted by Abrikosov [Q, for type II superconductors, upon increasing 
the applied field beyond the lower critical field Hd a type II superconductor undergoes 
a transition from the Meissner phase to the Abrikosov vortex lattice phase. In the 
vortex lattice phase the magnetic field penetrates the material and induces elementary 
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Figure 2: Crossection of a vortex line in type II superconductors, illustrating the 
relevant length scales discussed in the text. 



topological defects in the order parameter ip{r), which in three dimensions are vortex 
lines. In the Abrikosov phase the long-range order persists with ijj{f) = \^jJ(r)\e'^'^^'^ . 
However, instead of the phase being constant everywhere, it "winds" by 27r around 
a vortex line. The amplitude of ilj{r) is nearly constant everywhere away from the 
vortex core but is substantially suppressed inside the core of size ^ or ^2, for fields 
applied perpendicular to the a6-plane or 2;-axis, respectively. The magnetic flux per 
vortex is the flux quantum, 0o = hc/2e. When the vortices are well separated this flux 
is confined by circulating screening currents within a London penetration length A of 
the core. For an isolated straight flux line running parallel to an axis of cylindrical 
symmetry of the material the field is of the form. 



B(r) 



(3) 



J e"''/^ for r > A 
I -ln(r/A) for^<r<A . 

For the most common situation of magnetic field applied along the vortex 
line has structure on two length scales, ^ and A, and its crossection in a6-plane is 
illustrated in Fig. || 

In the Meissner phase and in the Abrikosov phase far away from the core the 
amplitude of the order parameter is a constant \ip{r)\ = that minimizes the 
pairing energy. Only the phase 0(r) varies substantially. Using this Ansatz for the 
order parameter ip{r) = ^/ple*"^^^ inside the Ginzburg-Landau free energy we obtain 
the London free energy which governs the phase and magnetic field variations. 



Pi 
2 



' d(f){r} 
dz 




27T 



+ ^ Vx0(r") - —A{f 

+ _L _ V X A{r, 
where and pi are the normalized eigenvalues of the superfiuid density tensor. 



Pi 



2mj 
= A„ 



Ps 



IGTr^A? 



(4) 



(5) 



where index i labels x, z and A^^ = Aj^ = A. 

The ratio k = X/^ determines whether the superconductor is type I (for which the 
Abrikosov phase does not exist) or type II according to whether k < 1/V2 or k > 
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l/-\/2, respectively. The high-Tc Cu02 materials are strongly type II, characterized 
by very large k At low temperatures the typical values of ^ < 15A and A > lOOOA 
leading to k > 100. 

The excess Gibb's free energy per unit of length for putting in straight, nonin- 



teracting vortex/flux lines along the z-direction is given by |30 



UIK — 



(6) 



A-nX J Air 

where H is the z-component of applied field H. This Gibb's energy vanishes at the 



lower critical field H = Hd [31 1, 



In mean-field theory both the coherence and the penetration length diverge as is 
approached from below, ^(T) ~ A(T) ~ (Tc — T)~^/^. Since this leads to k, that 
has a smooth behavior near we observe from Eq.^ that Hd vanishes linearly as T 
approaches Tc. 

For Hci < H < Hc2 vortices parallel to H begin to enter the sample with increasing 
applied field. The equilibrium number of lines is determined by the balance between 
the noninteracting contribution to Gibb's free energy, Eq.^, which is linear in A^, and 
the contribution nonlinear in A^, coming from the interaction between flux lines. The 
vortex-vortex interaction between two parallel straight lines separated by distance 
r > ^ is repulsive and is given by 

= ^ [Ko{r/X) - Ko{r/0] , (8) 
where Kq{x) is the modified Bessel function with the asymptotics, 

J^o[X) ~ j for x^O . ^ ' 

In Eq.p we have introduced a short distance cutoff for r < C,, with the interaction 
energy saturating at the value of excess energy for creating a doubly quantized vortex. 

In the absence of thermal fluctuations the flux lines crystallize into a hexagonal 
Abrikosov vortex lattice, a 3-d crystal of inflnitely long parallel vortex lines. The 
vortex lattice breaks continuous translational symmetry perpendicular to H and ro- 
tational symmetry about H, as well as the global U{1) symmetry. The flux line 
density n is directly related to the average magnetic fleld n ~ 1/a^ = {B{r))/(f)Q. 
Several values of lattice constants a, and the corresponding magnetic fleld strengths 
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Figure 3: Bitter pattern showing flux lines as they emerge from the surface of YBCO, 
at T = 4.2K. The applied magnetic field ifj, = 20 Gauss corresponds to 1/xm sep- 
aration between vortices. The very straight, somewhat smeared lines cutting across 
the lower two corners of the image are twin planes in the YBCO crystal. These twins 
locally orient the Abrikosov flux lattice. However, since the twins run at right angles 
to one another, the two resulting orientations of the hexagonal flux lattice are incom- 
patible. This results in the grain boundary in the flux lattice which runs from the 
lower center of the picture toward the upper left corner. 

are displayed in Table ^ For weak magnetic flelds resulting in flux line separation 
a ^ A, ^ the amplitude of the order parameter is y^p^ and 5 = away from the 

vortex cores. However as H ^ i/^, the amplitude l^/'l is suppressed while B grows 
continuous between vortices, until at Hc2 the vortex cores overlap and in the absence 
of fluctuations a mean-fleld transition occurs from the Abrikosov superconducting 
phase to the normal phase. 

Table 1: Relation between typical values of B and lattice constant a 

~B 20 Gauss 2 kGauss 20 Tesla 

a 1 /im lOOOA lOOA 



The existence of the Abrikosov vortex lattice for conventional superconductors has 
been demonstrated experimentally many years ago through various techniques that 
have now been applied to study vortex states in high-Tc superconductors 33|. The 



lattice can be imaged in real space through the Bitter decoration technique in which 
flne magnetic dust shows the location of flux lines as they emerge from the surface 
of the superconductor. Magnetic neutron scattering off B{r) due to the flux line 
lattice is another technique that has been successfully used. Figure |^ clearly shows 
an image of flux pattern on the surface of YBCO at T = 4:.2K. A 20 Gauss fleld 
corresponding to 1 vortex//im^ was used to obtain this Bitter pattern |3^. Figure^ 
shows the reciprocal-space image of the flux line lattice in NbSe2, a strongly type-II 
layered low-T^ ~ 7K superconductor. The image was obtained via small angle ~ 1° 
neutron scattering with the applied fleld of 8 KGauss at 5.2K. ||3^ 

Also, recently, new revolutionary electron holographic techniques have been used 
to image the motion of flux lines in thin Pb fllms. This probe, which can image flux 
lines in real time, can be used to study mixed states also in high-Tc superconductors, 
and can provide valuable information about the statics and dynamics of vortex lines 
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Figure 4: Reciprocal space image of the Abrikosov flux lattice, obtained using neutron 
small angle (~ 1") scattering from NbSe2, a strongly type II low-Tc superconductor 
(Tc ~ 7K). The image was obtained with 8 KGauss magnetic field applied along the 
c-axis, at T = 5.2K. The six hexagonally placed small spots are the magnetic Bragg 
scattering from the flux lattice. The large spot at the center is a nonmagnetic back- 
ground that includes small angle scattering from imperfections in the sample. The 
scattering at the higher-order Bragg spots is too weak to see under these conditions. 



Because we are interested in high-T^ superconductors we will consider strongly 
type II materials, A ^ ^, which means that A{r) is much "stiffer" than ip{r), corre- 
sponding to low carrier density or large effective mass. For H ^ Hd the 5-field of a 
vortex overlaps many other vortices, and one can often make a uniform field approxi- 
mation, B{r) = B =constant, and concentrate on the dynamics and variations of iplf) 
only. This regime can be formally obtained by taking limits e — > (0o oo) and 
H ^ oo with fixed vortex density n = {B)/(j)Q. This limit corresponds to uncharged 
rotating superfiuid ^He. 



2.2 Mean-field dynamics 



The mean-field dynamics within the infinite penetration length approximation is as- 
sumed to be accurately described by the time-dependent Ginzburg-Landau theory. 



dip{r,t) i2eV{r,t) 



5F, 



GL 



(10) 



dt h '^^ ^ 6ilj*{f,t) ' 

where for simplicity we have left out the Hall effect contribution. The first term in 
(|TUp is the Josephson term which leads to "rotation" of ip{r,t) in the complex plane 
as the voltage gradient £{r,t) = VV{r,t) accelerates the supercurrent js{r,t) = 
5FGL/5A{f,t) 



Re 



m Tic 



p.(V0-|l) 



(11) 
(12) 



We observe that the supercurrent is proportional to the net phase gradient. In the 
spirit of the two-fluid model of superfiuids the normal current is simply. 



(13) 



where m and (T„ must be replaced by their tensor equivalents for transport in an 
arbitrary direction, reflecting the intrinsic material anisotropy p7| . 
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Figure 5: Illustration of the origin of the Magnus force, that the transport current j 
(shown flowing in the xy-plane) exerts on the vortex line emerging from the plane. 

In the Meissner phase there are no vortices present, and therefore the phase is a 
true linear superconductor below critical current density jc- The transport current 
is however confined within a penetration length of the surface, js{r) ~ e~^/^, which 
limits the potential applications. For js > jc the resulting surface B field exceeds the 
local Hci, vortices begin to enter the sample. The ensuing vortex motion (in response 
the current-induced forces) generates dissipation and breaks down superconductivity 
as we explain in more detail below for the Abrikosov phase. 

To understand transport properties of the superconductor in the presence of vor- 
tices we must first examine the interaction between the supercurrent and the vortex 



line. The nature of this interaction is still a topic of much current research |38|. It 
is generally believed that there are two sources of interaction. One is the Lorentz 
interaction of the supercurrent js with the magnetic flux carried by the vortex line, 
resulting in the force fi = ^t{z) x j, where f{z) is a tangent unit vector to the 
vortex line. Since we are working in the large magnetic penetration length regime of 
nearly uniform B, the Lorentz interaction should not be very important. A more sig- 
nificant contribution in this regime is the interaction of the transport current with the 
screening supercurrents around the vortex. As is illustrated in Fig. |^, for the counter- 
clockwise supercurrent, currents add below and subtract above the vortex resulting 
in higher and lower kinetic energy, respectively. The vortex therefore experiences a 
Magnus force (analogous to the force that gives an airplane its lift) Jm = x j 

which moves it perpendicular to j, thereby reducing the kinetic energy and the phase 
difference. 

The evolution of the net phase difference between two points can be obtained from 
the dynamic Ginzburg-Landau equation ( p!0D by assuming that the amplitude of ipi^) 
is constant, 

= — 27r X (net rate at which vortices pass between the points) . (14) 

To maintain steady state therefore requires an applied voltage difference, 

AV ~ ^ ^ (^^^ ^^^^ which vortices pass by) , (15) 

which is an expression of the Josephson effect. 

The vortex motion in response to the uniform transport current is viscously 
damped by interaction of normal electrons in the vortex core with the ions, which in 
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Figure 6: Illustration of vortex lattice distortion due to the presence of non- uniform 
currents. As described in the text, since the total force and torque vanish, the lattice 
distorts but does not move as a whole. This therefore illustrates the vanishing of 
linear resistivity at finite wavevectors. 



turn results in dissipation. To maintain d.c. current an electric field £ needs to be 
applied proportional to the components of the transport current perpendicular to the 
vortices, 

^- = -F]—Pnj± , (16) 

-nc2 



where above proportionality constant is the Bardeen-Stephen flux flow resistivity p9 |. 

can be loosely interpreted as the ratio of normal to superconducting carriers, with 
only the normal ones contributing the normal-state resistivity p„, in the spirit of the 
two-fluid model of superfluids. 

Therefore for uniform current the transport is Ohmic perpendicular to vortices 
and is superconducting parallel to perfectly straight vortices (no thermal fluctuations) 

mm- 

The situation is quite a bit different for a nonuniform current ||4^ ; a configuration 



that is more physical since contacts and sample geometry that produce only uniform 
current are difficult to achieve. Consider the linear response of a vortex lattice (an 
elastic solid in a viscous medium) to a steady-state d.c. currents, neglecting surface 
effects. The steady state motion is determined only by the total force and torque. If 
the current is nonuniform and if the currents are fully balanced, the total force and 
torque vanish. Then the lattice is simply distorted by the current but no net steady 
vortex motion takes place. An example is depicted in Fig. ^ 
We can resolve j{r) into Fourier modes. 



j{r)= J d're''-^ j{k) , (17) 

and ask about the electric field response due to such a plane wave of current, 

S,{k) = p,,{k)Uk) , (18) 

where the summation over spatial indices is implied. 

As we already mentioned for the current parallel to the vortices, 

4(^) = . (19) 

The current perpendicular to the vortices couples to "phonon" modes of the vortex 
lattice, imposing static elastic distortions but no steady motion for k ^ 0. Therefore 
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we obtain a highly nonlocal resistivity with S{r) determined by j(r*) over the entire 
sample. 

pU^) - ^i^) ■ (20) 
We therefore conclude that in the absence of fluctuations, the vortex lattice super- 

— * 

conducts except ior k — and currents perpendicular to the vortices. Since in the 
normal state the resistivity is fully local, dissipating at all wavevectors k, Pij{k) — pn, 
in mean-field theory pij is discontinuous at the NS transition (at Hc2{T)). 



3 THERMAL FLUCTUATIONS 

We now examine how the mean-field results discussed in the previous sections are 
modified by thermal fiuctuations. In low temperature superconductors these fiuctua- 
tions are not very important except very close to phase boundaries and in thin films 
or wires. The small size of the critical region in bulk low-T^, materials is associated 
with the fact that the thermal length 

A, = . , (21) 

is much larger than any other length in the problem except very close to Tg where 
other lengths diverge. 

For the high-Tc superconductors simple considerations also allow us to construct 
the "Ginzburg criterion" , which determines when the effects of thermal fluctuations 
become important. This occurs when the ordering condensation energy in a coher- 
ence volume is equal to the thermal disordering energy, ~ /cbT. Stated 
equivalently in terms of material parameters, 

X^^l^Kr^^-^. (22) 

We therefore observe that for the highly anisotropic high-Tc superconductors the 
effective thermal length is reduced by more than 4 orders of magnitude by the com- 
bination of high Tc, large anisotropy ~ 10^ and large n ^ 10^. This therefore 
greatly expands the critical region in which mean-fleld theory no longer applies and 
fluctuation-corrected theory must be used. 

We first consider fiuctuations in the ideal case of disorder-free superconductors 
and then go on to study the effects of disorder at finite temperatures. 



3.1 Disorder- free samples 

The finite temperature effects can be incorporated into the statics by summing over all 
the configurations of the system weighted by Boltzmann weight 6"^°^/*^^^, instead of 
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just simply minimizing the Ginzburg-Landau free energy A], Eq.|l|. The effects 

of thermal fluctuations can equivalently be incorporated in the dynamic Ginzburg- 
Landau equation ( [TI1| ) by the addition of white noise ({r,t) of intensity proportional 
to ksT, as set by the fluctuation-dissipation theorem. 

Near the phase boundaries the fluctuations are most effective in disrupting the 
ordered phases, because there the stiffness of the order parameter vanishes. For 
example, for the Meissner phase p^^{T) ~ ^j^mf _ rp^^ Similarly, for the vortex 
lattice phase the shear modulus plays the role of the stiffness, which vanishes near 
the H^^ as 

{HM^_-Hy 
X 

The idea that the flux line lattice in high-Tc materials might melt into a liquid {fi 



l^iH) ~ . (23) 



0) of freely moving vortex lines was flrst suggested by Nelson [10|. Recently there has 



been much experimental evidence that in fact a large portion of the Abrikosov lattice is 
replace by the vortex liquid, all the way up to the mean-field H^2^ boundary where the 
superconductor crosses over to normal behavior |TT[ (see figures |I|(b), |^). Although 
a full theory of the melting transition does not exist, an estimate of the location 
of the melting boundary can be made using the Lindemann criterion. Much work 
has been directed to accurately describe the melting boundary, carefully taking into 
account the nonlocality of the elastic moduli of the vortex lattice However 
to avoid drowning in technical details we illustrate the idea ignoring the wavevector 
dependence of the elastic moduli. 

As is true for real ionic lattices, the phonon modes of the vortex lattice are ther- 
mally excited with mean elastic energy ksT /2. The mean square displacement of a 
vortex is a rough measure of the disruption of crystalline order by thermal fluctua- 
tions, 

{|'S(r,()p) = ld'k{\Mk,tW) , (24) 
'^-""''•-"^'^M: , (25) 

/ifc^ 



(26) 



As expected the displacement increases with temperature and as H^2^ is approached. 
Melting occurs when the displacement becomes some signiflcant fraction of the lattice 
spacing a, 

{\u{f,t)\Y^^CLa, (27) 
where cl ~ 0.1 is the empirically determined Lindemann ratio. 
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Figure 7: Some of the various newly proposed vortex phases are illustrated. The 
quasi-2d-supersolid vortex lattice phase is expected to be a distinct phase from 
the conventional Abrikosov 3d-lattice, and will exist at high 5-fields, where the in- 
ter-planar coupling is weak. 



The resulting vortex liquid has local short-scale superconductivity, but the highly 
mobile disordered vortices do not allow global phase coherence to set in. In this sense 
the vortex liquid is the same disordered thermodynamic phase as the normal state, 
exhibiting only the quantitative difference of having much larger fluctuations of the 
superconducting order parameter. In both of these normal-state regimes {ip) = 0, 
unlike in the truly superconducting Abrikosov lattice and Meissner phases. At H^2^ 
a gradual crossover occurs as well- defined vortices form, and conductivity and dia- 
magnetism increase due to the local superconductivity. However, no thermodynamic 
singularity occurs since no symmetry is broken as is crossed. 

The interaction between flux lines becomes very weak and falls off exponentially for 
separations r > A. Recently it has been argued that because of the weak interactions, 
at weak fields a sliver of vortex liquid will also intrude right above Hd fT^, |^, as 
is illustrated in Fig. 0(b). However the verification of the existence of the flux-line 
liquid in this portion of phase diagram has up to now alluded experimentalists. 

A variety of liquid phases have been suggested such as the entangled and disentan- 
gled vortex liquid |]10[ as well as the intermediate hexatic phase |jl2|, |l^ characterized 
by the quasi-long-range orientational order and short-range translational order in 
analogy with the 2(i-melting theory |T^. Recently, there have also been suggestions 
that in layered superconductors there is a possibility of two thermodynamically dis- 
tict vortex lattice phases. One is the conventional 3d Abrikosov lattice with vanishing 
linear resistivity. The second is a quasi-2(i vortex lattice phase, with finite linear re- 
sistivity, also known as the supersolid phase @, |^. In Fig. |^ we have presented a 
schematic phase diagram showing the expected location of these novel vortex phases. 
In the {B{f)) versus T phase diagram of Fig. 0, the Meissner phase is collapsed down 
to the T-axis. The nature of the novel hexatic and entangled liquid phases has been 
extensively discussed |1^, |13[. We now briefly describe the less known putative 



novel quasi-2(i lattice phase. 

In a usual molecular crystal an interstitial is a point defect that costs a finite 
amount of energy and in equilibrium exists at densities ~ ^-E,/kBT ^ ^^iq other 
hand in a three-dimensional vortex lattice a vortex cannot end inside a superconduc- 
tor (except on a magnetic monopole, which is highly unlikely) so an interstitial is a 
line defect with a constant free energy per unit of length. Its free energy is therefore 
proportional to the thickness of the sample. Since in the thermodynamic limit this 
defect energy diverges, at equilibrium we expect no interstitials in the vortex lattice 
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and magnetic flux per unit cell is exactly 0o- This situation is to be contrasted with 
thin film superconductors where the interstitials and vacancies are point defects, of 
finite energy and therefore are present in finite density at finite temperature. Con- 
comitantly the flux per unit cell does not equal to (po, and hence the number of flux 
points per unit cell fluctuates as the interstitial and vacancies move around. 

The crossover between the 2d and 3d behavior in layered Cu02 superconductors 



is described by the Lawrence-Doniach model [43| with free energy 



F = Y. {^2d [mm - j d^rRe [Jr^{m+i{^]] , (28) 

where we have taken B to be perpendicular to layers. The first term is the free energy 
of independent 2c? layers labeled by index i, and the second term gives local coupling 
between adjacent layers. 

At J = 0, the 2d uncoupled layers allow for vacancies and the interstitials with 
entropy dominating over energy. As the coupling is increased at some critical value 
Jc = J{Tc) the interstitial free energy per unit of length (tension) becomes positive, 
the system makes a transition to a 3(i-dominated behavior and the net interstitial 
density vanishes. Since the effective interlayer coupling is actually J = Jo? ~ J / (B) 
the system becomes more 2c?-like at higher fields. At large fields (small J) the inter- 
stitials form a highly flexible and mobile line-liquid existing inside a 3d vortex lattice 
with a true crystalline order. Since the interstitials are free to move, this supersolid 
is not superconducting for the reasons explained in previous sections. The resulting 
phase diagram is illustrated in Fig. ^ where we have plotted (B) versus T. 

3.2 Nonlocal resistivity at finite temperature 

As we already discussed in Sec. 2. 2, the transverse flux-flow resistivity of a vortex 
lattice is highly nonlocal, 

S{f) = J d^r'p{f-f')j{f^) , (29) 

and in fact in the absence of plastic flow (possible in very clean samples) we expect 
Piatticei^ ~ 1/volume. This form is to be contrasted with the completely local resis- 
tivity in the normal state, Pnormai ~ 6^^\r). In the presence of thermal fluctuations 
the flux-line lattice melts over a large portion of the HT phase diagram. The resulting 
vortex-line liquid consists of flexible, highly mobile flux lines that can easily respond 
to the transport current and therefore generate finite flux-flow resistivity. Unlike the 
vortex lattice, in the vortex/flux-line liquid phase there is no perfect correlation in 
the motion between different, separated vortex lines. As the lines traverse the sample 
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Figure 8: Schematic of the experimental setup to measure nonlocal resistivity of the 
vortex liquid phase. In contrast to the normal state, transport in the vortex liquid 
will generate Vtop = Vbottom- 



along the applied magnetic field (say in the 2;-direction) in a distance z they will 
wander transversely a distance {zksT /ty/"^ , where the line tension e = 0o-f^ci/4vr. 
Even in thin samples a typical line will deviate away from its average direction much 
further than the average interline separation set by the applied magnetic field (see 
Table |l]). Therefore, the flux lines experience many collision with their neighbors 
entangling with many distant ones, much like directed polymers in a solution. Al- 
though no rigorous calculation of the fiux-line crossing barrier energy is available, 
simple estimates give ~ 50k bT [|1^, far away from Hc2, where it is expected to 
vanish as ~ \n[H^2^ / H]/ \ia[H^^ / Hd] [0]. We expect that these frequent flux-line 
interactions and entanglements will lead to a large transverse viscosity of the vortex- 
line liquid |]rB| which replaces the perfect correlations controlled by the finite shear 
modulus in the Abrikosov lattice. Besides these transverse viscous correlations, the 
connectivity of the vortex lines leads to an even stronger nonlocality in the viscosity 
of the line-liquid along the direction of the applied field. Since the viscosity is a 
measure of the response of velocity to an applied force, it is directly related to the 
fiux-fiow resistivity. The fiux-fiow resistivity of vortex liquid should therefore exhibit 
a nonlocal behavior in-between that of a completely local one of the normal state 
and a completely nonlocal one of the vortex lattice. As the superconductor is cooled 
into the vortex liquid, p{r) becomes more and more extended, especially parallel to 
vortices. 

The nonlocality of the resistivity in the liquid phase has been recently studied 
experimentally W\ and theoretically |^7[ in the geometry depicted in Fig. |^. 

The transport current at the top surface of the superconductor exerts a transverse 
force on top ends of vortex lines. Since the vortex lines are well connected and rarely 
break, the current at the top of the sample pulls on the whole line. This results in the 
vortex motion and therefore electric field far away from the region where the current 
actually flows (top). Since the voltmeter measures the average rate at which vortex 
lines pass between the two contacts, excluding the possibility of line breaking we have 
Vtop = Vbottom, with the vortex velocity proportional to the total current integrated 
from the top to the bottom of the sample. This situation is to be contrasted with 
the normal state local resistivity where the Vtop > Vbottom because current density is 
higher at the top. In experiments of Safar, et al. [K] on 35fim thick YiBa2Cu30j 
crystal, in a regime where vortices have a good integrity over a distances larger than 
the sample thickness, the effect of nonlocal resistivity has been observed. 

Supposing that the Fourier transform of the conductivity is an analytical function 
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of A; = {k±,kz) for small k, nonlocal resistivity of the vortex liquid can be 
phenomenologically described as [l^, |T6| 



^ ^ It 12^ ,2 ' ^30) 

7 + Vab\k±\^ + Vzki 

where 7 is the local friction (related to the local conductivity of Bardeen and Stephen), 
and rjab, tJz are the vortex liquid viscosities leading to nonlocality in the flux-flow 
resistivity. This model can be used to solve for current and voltage patterns in a 
realistic experimental geometry by solving 4th order partial differential equations. 

3.3 Static disorder 

In real high-temperature superconductors disorder plays an essential role. As we have 
already briefly described, it modifies the HT phase diagram, but more importantly, 
the disorder tends to pin vortices in place, drastically reducing the fiux-fiow resistivity, 
as we now describe in more detail. 

The most important contribution of static disorder leads to a spatially varying 
transition temperature, Tc(r). The disorder effects can therefore be incorporated into 
the Ginzburg-Landau description by allowing the a parameter to be a function of 
position, 

aif,T)^a{r)(T-T,''^{r)) , (31) 

where the properties of the function depend on the nature of disorder. Chemical 
and structural imperfections in the material naturally lead to random a(r^ with only 
short-range correlations. Correlated random static disorder is also possible and might 
arise from grain or twin boundaries, screw dislocations, artificially created ion colum- 
nar tracks [^, 0, or from an epitaxial multilayers and e-beam written patterns. 
Finally the crystal structure itself, such as for example the CUO2 layers of the high- 
Tc superconductors, provides a natural modulation of Tc(r) describable by a periodic 
a[r). 

Since the core of a vortex line is in a normal state, the energy is minimized when 
the vortex is located in the region of weakest superconductivity, where a{r) is greatest 
and minimum amount of condensation energy is lost. This leads to pinning of a vortex 
line by normal impurities, dislocations, grain/twin boundaries and to the intrinsic 
pinning in between the CUO2 planes (where the superconductivity is weakest). 

The phase of the wavefunction cannot be measured and therefore the disorder 
cannot directly couple to the phase 0(r). We therefore conclude that the Meissner 
phase is stable to disorder, except near its phase boundaries. 

Both the strong and weak pinning limits have been previously studied. In the 
strong pinning limit the pinning energy exceeds the vortex-vortex interaction energy 
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and leads to a fully disrupted vortex lattice. The Ginzburg-Landau free energy is 
minimized by a random complex function ip(r) = ipQ{r)e^'^^'^\ 

It has been known for some time that the Abrikosov phase is unstable to intro- 
duction of disorder ||^ and the long-range translational order is destroyed beyond a 
disorder-determined Larkin-Ovchinnikov length L^q- K has been recently proposed 
that the resulting new low temperatures and fields phase possesses long-range order 
with a sample and disorder-specific random vortex pattern. For the case of point (un- 
correlated) disorder the equilibrium phase is the vortex glass with the name chosen 
by the analogy with spin glass, where the ground state order parameter is random but 
frozen with long-time correlations |51|, For correlated disorder such as columnar 



pins, an anisotropic version of the vortex glass has been proposed and dubbed Bose 
glass 1^ by the isomorphism of the resulting theory with bosons on random substrate 
Here we will confine our review to the isotropic vortex glass, with the analysis 
easily extendible to Bose glass. 

The stability of the vortex glass phase to thermal fluctuations depends on the 
nature of the low- lying states. By analogy with the scaling theory of spin glasses |l53| , 
the excitations are "bubbles" inside a ground state and correspond to the displacement 
of one or more vortices a distance L in a region of linear size L. The average energy 
of the low-lying excitation is expected to scale with a power of L 

~ TL^ . (32) 

Since for a circular loop, 9i = 1, characteristic of excitations in the Meissner phase, 
we expect that 6 <1. 

The stability of the vortex glass phase clearly depends on the sign of the 9 ex- 
ponent. If ^ < then the large-scale excitations cost very little energy and will be 
present at T > at equilibrium. This turns out to be the case in two dimensions 
{d = 2). The correlation length then scales as a power law with temperature 

ivG ~ T-VI^I , (33) 

with the long-range order only present strictly at T = 0. 

In three dimensions no satisfactory analytical analysis is available to date, pri- 
marily because of the difficulty associated with treating the dislocations (although 
some arguments have proposed that the dislocations might be irrelevant). Numerical 
studies of highly simplified models treat the 3d vortex glass phase in the limit A ^ oo 



and find < 6' < 0.3 |^3|, |2^, These works point to the stability of the vortex 
glass phase in = 3, for T > 0. 

The stability of the vortex glass phase in three dimensions to thermal fiuctuations 
in the presence of screening (A finite) is still an open question. Since the interactions 
crucial for the formation of the vortex glass phase become exponentially weak for 
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r > A, it appears that screening tends to weaken the stabihty of the vortex glass. 
Experimentally this regime is only probed when ^vg > ^ which corresponds to resis- 
tivities well below those studied in experiments to date. 

In the case of weak disorder a topological glass phase might result, with properties 
intermediate between that of the vortex lattice and the vortex glass. As we have seen 
vortex glass is a phase highly correlated in time, but spatially it looks like a frozen 
vortex liquid with many dislocations present. In the vortex glass the strong disorder 
destroys the topological order and dislocations proliferate. On the other hand, in 
the absence of disorder, the vortex lattice has both positional and orientational long- 
range order and is therefore dislocation-free at large scales. It is possible that in 
the presence of weak disorder an intermediate topological glass forms in which the 
topological order of the vortex lattice survives. 

The resulting topological glass phase is then describable as a dislocation-free se- 
lastic medium (lattice) in a random potential. The displacement of a vortex from a 
position r = (x, z) due to disorder is described by u(r) ± z, with vortices parallel to 
z. The Hamiltonian in d dimensions is then given by 

H = j d'^r Vu(f)- E ■Vu(f) 

+ T.JdzV ({xi + u(xi, z)}, z) , (34) 

i 

where Xj is the d—1 dimensional position of the i-ih line in the space transverse to z 
and the sum is over discrete vortex lines. The first term, above, is the elastic energy, 
quadratic in the displacements u(r), described in terms of the elasticity tensor E 
[0, |T^. The second term is the contribution due to the pinning disorder, described 
by a random potential V{r) with spatial correlations that depend on the type of 
disorder. 

The ground state without dislocations has been studied analytically and is char- 
acterized by displacement correlations that are logarithmic in position 

[u(fi) - u(f2)]^ ~ In |fi - , (35) 

for large \fi — [Q, Weak disorder destroys long-range positional order re- 
placing the Bragg peaks by power-law singularities characteristic of the quasi-long- 
range order. Because of the weakness of the logarithmic growth of displacements, the 
dislocation density vanishes and the topological order remains, consistent with the 
assumption of the model. Full long-range orientational order, however, still survives 
in this randomly pinned vortex lattice. 

For the topological glass to exist as a distinct phase from the strongly disordered 
vortex glass it must be stable to dislocations. To establish this stability we look 
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Figure 9: Illustration of a domain wall of nonlinear distortion resulting from shifting 
a patch of topological vortex glass by a lattice vector a. The single valued vortex 
displacement Uj(r) inside the wall differ from the distortions outside the patch, Uo(r) 
by a. The displacement field jump of a is concentrated on the domain wall. 

Figure 10: An addition of a dislocation line together with a negative energy domain 
wall might lower the energy of the topological glass, depending on the value of the 9^ 
exponent. This would result in distablization of the topological glass toward a fully 
disorder vortex glass. 

at the low-lying excited states that disrupt the order of the topological glass. An 
appropriate low-energy excitation is a patch of size L moved over by one lattice 
spacing a with respect to the rest of the vortex lattice. The elastic nonlinear distortion 
is Uo(r), outside the patch and Uj(r) inside the patch, with the displacement jump 
|uj(r^a//) — ^oifwaii + an) | = a, concentrated on the wall bounding the patch (see Fig. 
1^). The elastic energy is unchanged in the interior and exterior of the patch. The 
average total energy is increased due to the elastic strain by a and for optimized wall 
position scales as 

E,{L) ~ TL^ . (36) 

Since the elastic energy is concentrated on the wall, the average total energy can at 
most scale as the area of the wall, which puts an upper bound on 

e^<2, (37) 

where the optimization of the wall position can lower 6*^, below 2. 

Although the total average energy of the wall is positive, large sections of the 
wall can have a negative energy that scales in the same way as the average energy in 
Eq.^. A section of the wall with negative energy can be created by an addition of a 
dislocation line, as illustrated in Fig. |TD|. Upon encircling the dislocation line Uj — uq 
must change by a Burger's vector b, a change that is concentrated on the domain 
wall. Therefore the energy of adding a dislocation equals to the local core energy plus 
the wall energy. On scale L the system might prefer to create a dislocation line with 
the negative energy domain wall compensating the core energy cost 

^dislocation ^ L — T ^ . (38) 

Hence we expect that the system can lower its energy by adding dislocations on 
large scales if 6'„, > 1 and the topological glass phase will always be unstable to the 
vortex glass without topological order. On the other hand if < < 1 the topo- 
logical glass phase should be stable against dislocations for weak disorder. To date 
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the precise range into which 6^ falls, and therefore the stability of the topologically 
ordered, dislocation-free vortex glass phase is an open question. 

4 NONLINEAR RESISTIVITY IN SUPERCON- 
DUCTING PHASES 

Vanishing linear resistivity is probably the most natural and practical distinction 
between a superconducting and normal phase. As we have already seen, the vortex 
liquid has nonlocal but nonvanishing linear resistivity. On the other hand a vortex 
lattice can be easily pinned with few strong inhomogeneities or even naturally pinned 
by the boundaries. In this pinned state we have seen that the vortex lattice has a 
true superconducting response to a uniform current. 

Conventional theory of transport in disordered vortex states, based on the ideas 
of Anderson and Kim describes the resistivity in terms of thermally activated 
motion of independent vortex bundles over the finite energy barriers Vp, introduced 
by the disorder ||^, 0. The size of the bundles is assumed to be finite, on the or- 
der of Larkin-Ovchinnikov length Llq, beyond which the lattice order is destroyed 
by disorder 0]. The transport current introduces an average tilt to the random 
landscape potential and therefore an overall thermal drift of vortex bundles in the 
direction of the tilt. The vortex drift destroys the phase coherence of superconduct- 
ing wavefunction, resulting in finite vortex-flow linear resistivity of Arrhenius form, 
pi ~ exp{—Vp/kBT). Therefore this conventional theory maintains that at finite tem- 
peratures, the disordered low temperature vortex phase is not really superconducting, 
so is qualitatively identical to the non-superconducting vortex liquid regime. 

The difference between the conventional picture and the transport in the vor- 
tex glass phase comes from the strong inter-vortex interactions not included in the 
Anderson-Kim picture. Upon cooling a sharp thermodynamic transition takes place 
from a non-superconducting vortex liquid to a distict superconducting vortex glass 
phase [|l^. The vortex motion in the vortex glass phase is still described by the 
thermally activated motion, except that the strong correlations lead to bundle and 
therefore barrier sizes that diverge as j 0. We now examine the vortex glass 
phase along with Meissner phase, and show that these phases possess only nonlinear 
resistivity and are therefore true linear superconductors. 

In the Meissner phase it is the nucleation and growth of vortex loops, and in the 
vortex glasses phases it is the collective motion of already present vortices in response 
to a transport current that are the excitations that lead to energy dissipation. Since 
the motion of a bundle of vortices can be described by the creation and growth of 
vortex loops, in both phases the low-lying excitations are vortex loops superimposed 
onto the ground state background of each phase. We can use the scaling theory of 
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these low-lying excitations to analyze the resistivity both in the flux-free Meissner 
phase and the frozen vortex glass phase. 

For the well understood Meissner phase we know the nature of the excited states 
(vortex loops) and therefore can predict the exponents in addition to the scaling form 
of resistivity. Unfortunately, for the vortex glass neither the ground state nor the 
low-lying excited states are well understood. This lack of understanding however, 
can be neatly packaged into a set of exponents that enter the scaling theory of the 
nonlinear resistivity. Some attempts have been made to apply a collective pinning 
theory to calculate these exponents. However, these theories ignore the existence 
of dislocations and therefore breakdown of long-range translational order on length 
scales L > Llo [i^- It is likely that dislocations will have an affect on the value of 
these exponents. 

As was described in the previous section, the energy of low-lying excitations scales 
as TL^ , with 6 being specific to the type and shape of the excitations and the phase- 
background in which they are created. For example for the isotropic Meissner phase 
the lowest energy excitations are circular loops characterized by superfluid density 
stiffness T ~ and 6 = 1. The transport current-generated Magnus force couples to 
the area A projected onto the plane perpendicular to the current 

A{L) ~ AqL^ . (39) 

The exponent r has a lower bound of 2, saturated by a circular loop in the Meissner 
phase with Aq = 1 and r = 2. The energy supplied by the Magnus force is 

E,{L) ^fMA{L), (40) 
~JL^ (41) 

where j is the transport current density. Balancing the elastic energy of the excitation 
with the current energy we obtain the average size Lj of the lowest energy excitation 

L,~r^/^^-^ (42) 

that are produced when t > 9. 

In order to reach these low-lying states the system might have to overcome barriers 
by moving through less favorable regions. In general the size of the energy barriers 
will scale as 

Eb = AL^ , (43) 

where for loops in the Meissner phase ip = 6 = 1. In general, however, the barriers 
are usually larger than the final energy state, and therefore 6' is a lower bound for ip, 



i.e. ip > 9. These two energy scales are schematically depicted in Fig.ITT 



23 



Figure 11: Local random potential due to impurity disorder is characterized by the 
typical size of its barriers, scaling as L^, and by the energy difference between ground 
state and next low-lying state, scaling as . 



The rate of energy dissipation via the production of the low-lying excitations is 
governed by the Arrhenius law ~ exp —AL^/ksT . The resulting nonlinear resistiv- 
ity due to the vortex motion over the barriers is 

where ^ = ip/{T — 6). The excitations with the smallest /x dominate with the upper 
bound again set by the vortex loops as in the Meissner phase where f^Meissner = 1- 
We immediately observe that for t > 6 the barriers diverge as j ^ and the linear 
dc resistivity vanishes 

fir"-"- <^^> 

Another consequence of Eq. 44 is a slow decay of non-equilibrium screening super- 
currents that can be detected by measuring the relaxation of the associated magne- 
tization M(t), 

dMjt) dj{t) 

^ oc ^ oc £ij) , (46) 

where S{j) is the field needed to prevent j{t) from decaying. Combining ( ^6[) with 
( PI) leads to a slow logarithmic decay at long times, 

j(t)^Jo[ln(tAo)]"'/" , (47) 

with to as the microscopic time of order 10^^ — 10~^^ seconds. 

Recent transport experiments on YBa2Cu30i thick films (in 3d regime) find non- 
linear resistivity behavior described above with fi ^ 0.2 — 0.3 which sets in when the 
temperature is lowered past a well defined transition temperature Tyg (see next sec- 
tion) |19, 20]. Later much more sensitive magnetization decay experiments on YBCO 



films have found a scaling behavior in agreement with Eq.^, with ^ 1/3 |^, ^ 
These experiments therefore provide confirmation of the vortex glass picture. 



5 PHASE TRANSITIONS 

Having discussed the nature of the phases that occur as a result of thermal fluctuations 
and disorder we now turn our discussion to the phase transitions between these new 
phases. 
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5.1 First-order vortex lattice melting transition in pure crys- 
tals 



In the absence of disorder the low temperature and low field phase is the vortex 
lattice, which melts into the vortex liquid upon increasing field and/or temperature. 
This transition has been extensively studied experimentally, numerically and analyt- 
ically. Transport experiments on clean YiBa2Cu30'j crystals show an abrupt (with 
milliKelvin resolution), hysteretic resistance drop to zero at the melting transition 
[|60[] . This experimental work provides solid evidence for the vortex lattice melting 
transition being first-order. The melting transition has not yet been detected thermo- 
dynamically (i.e. by measuring latent heat or magnetization jumps) due to smallness 
of the vortex latent heat, estimated to be < ksTM per vortex per CuO layer. The 
vortices are very dilute and the specific heat due to electrons and phonons of the 
underlying solid are quite large at Tm ~ 70K in comparison and dominate over the 
vortex lines' contribution. Also since the high-T^ materials are highly anisotropic and 
therefore are in the quasi-2d regime, the first-order transition is expected to be weak. 

Evidence from the simulations of London theory (amplitude oiijj{r) is a fixed con- 
stant) on a lattice also clearly points toward the first-order transition [^. Further 
evidence is provided by the analytical work of Brezin, et al. [^|. There an e = Q—d ex- 
pansion is carried out which shows that the mean-field theory second-order transition 
at is destabilized by fiuctuations. The resulting runaway of renormalization- 

group fiow is taken as an indication of fiuctuation driven first-order transition |^ 



Finally based on general symmetry grounds Landau theory for melting has a term 
cubic in the order parameter and therefore leads to first-order melting transition in 
agreement with usual melting of 3d crystals [pi]. 



5.2 Stability of first-order transition to disorder 

We have previously described the strong infiuence disorder has on vortex phases, for 
example converting the vortex lattice to a vortex glass. We now examine the stability 
of the first-order melting transition to weak disorder using Imry-Wortis general scaling 
arguments [|65|| . 

Near the melting transition, disorder's local and random preference for the solid 
(vortex glass) or vortex liquid phase can be described via Landau theory with a 
random local variations in Tm (see Eq.^). An increase in free energy due to a creation 
of an island of one phase inside the other will have a positive interface contribution 
and a favorable negative bulk contribution due to the disorder 

5F ~ aL'^-^ - AL'^/^ , (48) 

where a is the tension of the c?— 1-dimensional interface and A is the disorder strength. 
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Figure 12: The nature of phase transitions at weak (l**-order) and strong (2"''^-order) 
magnetic fields is illustrated. The diagram is based on combination of theoretical 
arguments and recent experiments results (see text). Tricritical point separating two 
kinds of transitions at approximately 10 Tesla (depending on the strength of disorder) 
is indicated. 



The negative disorder contribution is a sum of random energies (with zero mean) and 
therefore by central limit theorem scales as a square root of number of impurities 
in the region of size L. For d > 2 and weak disorder (small A), the free energy of 
the island is positive and their density will be strongly suppressed. Therefore, there 
will be a stable two-phase coexistence with the phase transition remaining first-order. 
On the other hand for d < 2 the disorder contribution will dominate the interfacial 
energy and system will create interpenetrating islands of one phase inside the other. 
These arguments therefore suggest that the melting transition of the vortex glass into 
a vortex liquid will become a continuous transition as the effective dimensionality 
decreases and the strength of disorder increases. 

Recent experiments on untwinned, single crystal Y Ba2Cu^0'j by Safar and cowork- 
ers find a behavior that can be reconciled with the above theoretical picture. They 
find a hysteretic, abrupt, first-order transition at weak applied fields H which becomes 
a continuous transition at higher H. It is also observed that the first-order regime 
increases as pinning decreases with the tricritical point moving in a range 5 — 15 
Tesla. Because increasing H is equivalent to lowering Tm and making the system 
more two-dimensional, it is believed that increase in H corresponds to an increase 
in the effective strength of disorder. A schematic of experimentally observed phase 



diagram is illustrated in Fig. 12 



5.3 Scaling theory near continuous transition 

The transition from the superconducting Meissner phase to the non-superconducting 
vortex liquid is continuous since the density of vortex lines (related to the order pa- 
rameter) increases smoothly at iJ^i- Also we have argued in the previous sections that 
at high applied fields and/or strong disorder the transition from the superconducting 
vortex glass phase to the vortex liquid is second order. As a continuous transition 
is approached a correlation length ^(T), characterizing the length scale over which 
correlations in the order parameter exist, begins to grow, diverging exactly at the 
transition T = Tc with a critical exponent u 

aT) = ^o\{T-T,)/Tr (49) 
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To the correlation length corresponds a correlation time r over which the correlations 
of size ^ decay. In the critical region the correlation length is usually larger than any 
other physical length in the problem and r scales with ^ as r ~ Based on these 
observation a continuous transition can be studied using a scaling theory, where our 
ignorance about the transition can be packaged into a small set of unknown universal 
critical exponents. Above we have confined our analysis to an isotropic situation 
with a single diverging correlation length. For anisotropic transitions, such as for 
example the transition to anisotropic Bose glass (due to correlated disorder) our 
analysis can be easily generalized by keeping track of scaling with correlation lengths 
in all inequivalent directions . 

The scaling theory of nonlinear resistivity can be established by identifying the 
scaling of the characteristic current j and the electric field S with the correlation 
length ^, and constructing a dimensionally-correct combination from these physical 
quantities. By definition, in the free-energy a supercurrent couples to the gauge 
invariant superfluid velocity and therefore the free-energy of the superconductor will 
have an additional current contribution, 

...^-/.^.^J.(v,-^). ,50) 

The long-range correlations in the superconducting order parameter ijj are lost beyond 
the coherence length scale. Equivalently, the phase changes by ~ 27r on the scale 
^ and therefore |V0| ~ Svr/,^. From above equation (and dimensional analysis) the 
characteristic value of the vector potential A is 

1^1 ~ f , (51) 

which when combined with equation S = —c~^dA/dt gives the scaling of the charac- 
teristic electric field 



Similarly, from Eq.^ we find that the characteristic supercurrent density in the cor- 
relation volume is 

cksT 

In three- dimensions above equation has a simple interpretation as a requirement that 
in the critical region the Magnus energy, in which current j couples to the projected 
area be of the order ksT. 
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Combining above equations we obtain a scaling form of the nonlinear resistivity 
for uniform current, valid near the NS transition 

M-'-^e-'--F.(^,e-), (54) 

where the scaling functions -F+, F_ describe the nonlinearity in resistivity for T > 
and T < Tc, respectively. Above scaling equation can be recast in a more convenient 
form in terms of T and j using Eq.^ 

where, 

s =u{2 + z-d) , (56) 
C =^{d-l), (57) 

and G±{x) are some other scaling functions. 

The frequency dependent linear resistivity is finite everywhere and the scaling 
analysis similar to the above leads to 

PiiT,u)^^'-^-'G±iujC) ■ (58) 

The properties of the phases above and below gives us the limiting forms of 
these scaling functions. Above the material is normal and the behavior is Ohmic 
at small j. Hence, 

G+(x) constant, for x — > , (59) 

with the dissipation becoming nonlinear for j ^ jni ~ (T — Tc)'^. The linear resistivity 
for j <^ jni is therefore expected to vanish as the transition is approached 

Pi{T)^{T-T,y . (60) 

On the other hand for j <^ jni, i.e. in the limit of T — at j ^ 0, we expect 
nonlinear resistivity to be finite at T^. This requires 

G±{x) x"^^ , for X ^ oo , (61) 

so that the divergences in [T — T^) cancel out. At we find nonlinear current- voltage 
relation 

£(J)~J^+^/^ (62) 

that is a power-law just like in the Kosterlitz-Thouless transitions p7 



28 



Finally, below Tc the power law iv-characteristics is replaces by a stretched expo- 
nential, described in Sec. 4, Eq.(^), 

^.(x) ^ e-<^/^^ forx^O. (63) 

The scaling behavior emerging from above equations have been seen in dirty, heav- 



ily twinned Y Ba2Cu30f samples |]T9|, ^ with an impressive scaling over 4 decades for 
the linear resistivity above and over 2 decades for Pni{j) at the vortex glass-liquid 
NS transition. They find s ~ 6.5 and 2z/ 4, critical exponents that are consistent 



with the values obtained from numerical studies [p3| , p5 |. 

Just as for normal metals skin depth for magnetic fluctuations (photons) at fre- 
quency uj is 

A(^)~(^)'''. (64) 
In the superconducting phase pi{uj) ~ iuj/ps and therefore 

A(^ - 0) ~ ^ (65) 

Ps' 

Since near Tc we found that pi ~ and the characteristic frequency ~ we 
conclude that A ~ ^^^"^ and the Ginzburg-Landau parameter is 

= ^ ~ r'^' . (66) 

This means that even if the superconductor is characterized by type II "bare" prop- 
erties, as the transition is approached and C,(T) oo, the superconductor becomes 
effectively type I. This crossover was first investigated analytically by Halperin and 
coworkers who found that the charge is strongly relevant at the NS transition. Al- 



though initial indications based on renormalization group calculations suggested that 
the fluctuations drive the transition to be first-order, subsequent numerical analysis 
points to the transition being that of the inverted XY- model [^, although the story is 
not fully settled. Recent self-consistent analytical methods also suggest a continuous 
second-order transition, but with different exponents |7D|]. If the bare kmf = 100, the 



crossover to this new critical behavior will occur very close to the transition, when 
^(T) has grown by a factor of 10^ from what it was when the critical regime was 
first entered. Unfortunately, this is experimentally too close to the transition to be 
detectable by the presently available techniques. 
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6 CONCLUSION 



As we have seen, the vortex states of high-Tc superconductors are highly correlated 

and strongly fluctuating disordered systems that therefore share much in common 
with other systems well studied and of much interest in condensed matter physics. 
Many analytical tools have been taken over from widely different fields of polymer 
physics, magnetic spin-glass systems, strongly interacting electrons and bosons in ran- 
dom environment, to name a few. As we described in these lectures, much progress in 
understanding these vortex states has already been accomplished, and a theoretical 
picture consistent with experiments and numerical simulations has emerged. Because 
the physics of the novel high temperature superconductors incorporates the most de- 
veloped and exciting branches of modern condensed matter physics we expect many 
further exciting developments and surprises in our understanding of the phenomenol- 
ogy of high-Tc superconductors. 
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